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ANALYSIS OF LAMINATED, COMPOSITE, CIRCULAR CYLINDRICAL 
SHELLS WITH GENERAL BOUNDARY CONDITIONS 

By-S. Sr ini v as* 

Langley Research Center 

SUMMARY 

This report develops (1) a refined approximate theory for the static and dynamic 
analyses of finite, laminated, composite, circular cylindrical shells with general boundary 
conditions; (2) an exact three-dimensional analysis of simply supported, laminated, com- 
posite, circular cylindrical shells; and (3) athin-shell theory for laminated, composite, 
circular cylindrical shells. 

In the refined approximate theory the displacements are assumed to be piecewise 
linear across the thickness and the effects of transverse shear deformations and trans- 
verse normal stress are included. A variational approach is followed to obtain the 
governing differential equations and boundary conditions. A general solution of the govern- 
ing differential equations is also presented. The analysis of finite laminated shells with 
general boundary conditions involves satisfying the boundary conditions by making use 
of the appropriate part of the general solution. By using the refined approximate theory, 
the edge boundary conditions can be properly satisfied. 

In the exact three-dimensional analysis of simply supported, laminated shells, each 
ply is treated as a homogeneous cylinder. The displacements are chosen to vary trigo- 
nometrically in the axial and circumferential directions. The three governing partial 
differential equations are reduced to three coupled ordinary differential equations with 
the radial coordinate as the independent variable. The three coupled ordinary differential 
equations are then solved by use of the Frobenius method to obtain the variation of displace- 
ments in the radial direction. By satisfying the interface and exterior surface conditions, a 
set of simultaneous algebraic equations is obtained. For free vibration, the determinant of 
the coefficient matrix is equated to zero, and the solution of this characteristic equation 
yields the frequencies. The analysis and results are applicable to wave propagation in 
infinite shells, since the simple boundary conditions simulate the axial nodes in infinite 
shells. 


*NRC-NASA Resident Research Associate. 



The thin-shell theories for laminated, composite, circular cylindrical shells that 
are available in the literature contain some approximations in addition to the basic 
assumptions of thin-shell theory. The present thin-shell theory contains no additional 
assumptions. 

The results obtained from the refined approximate theory and thin-shell theory 
are compared with the exact results for simply supported, laminated, composite circular 
cylindrical shells. The refined approximate theory is found to be very accurate even 
for thick shells with short nodal distances. In contrast, the thin-shell theory is found 
to be reasonably accurate only for thin shells with moderate nodal distances. Frequency 
calculations using the refined approximate theory and thin-shell theory took considerably 
less computer time than those using the exact three-dimensional theory. 

INTRODUCTION 

Structural application of laminated composites in which each ply is made up of 
different materials is ever on the increase. As a consequence, analysis of such structures 
is gaining importance. The present paper deals with the static and dynamic analyses of 
laminated, circular cylindrical shells made up of composite materials. 

Analysis of laminated circular cylindrical shells has drawn considerable attention 
in recent times. If any of the plies is fiber reinforced or honeycomb, the normal proce- 
dure is to treat the ply as an equivalent homogeneous material having orthotropic or 
anisotropic properties. Dong (ref. 1), Bert, Baker, and Egle (ref. 2), and Stavsky and 
Loewy (ref. 3) have given differential equations for the dynamics of laminated cylindrical 
shells and solutions for free vibrations when the shell is simply supported. In these ref- 
erences the differential equations are based on different variations of thin-shell theory. 

The basic assumptions of thin-shell theory are (1) radial displacements are constant 
across the thickness, (2) axial and circumferential displacements vary linearly across 
the thickness and, (3) transverse shear deformations and transverse normal stress can 
be neglected. In references 1 to 3 some additional approximations are made. 

Nelson, Dong, and Kalra (ref. 4) have analyzed the special problem of free vibrations 
of simply supported circular cylinders by following the Ritz technique. Using the theory 
of three-dimensional elasticity, Armenakas (ref. 5) has analyzed the problem of wave 
propagation in two- layered circular cylindrical shells of infinite length and made of iso- 
tropic materials. For isotropic materials, the governing differential equations of elasticity 
can be easily solved in terms of displacement potentials; the variation of displacements 
in the radial direction, when determined, is in the form of Bessel functions. However, a 
similar approach is not possible if the material is orthotropic. 
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For laminated composite shells with general boundary conditions, the literature does 
not appear to offer any analysis which takes into account general variations of the displace- 
ments across the thickness. Such an analysis would automatically take into account the 
effects of transverse shear deformations and transverse normal stress. An analysis of 
this type is needed for the following reasons: 

(1) The edge boundary conditions can be properly satisfied only if the variation of 
cli gplacements across the thickness is unrestricted. The stresses and displace- 
ments close to the boundaries obtained from the thin- shell theory are highly 
inaccurate. 

(2) For thick shells, especially composite shells, the thicknesswise variation of 
displacements is more general than the restricted distribution assumed in the 
thin- shell theory. Therefore, a general variation of displacements is necessary 
for a better estimation of frequencies of free vibrations, stresses, and displace- 
ments even in regions far from the edges. 

An exact three-dimensional elastic analysis of laminated shells is computationally 
impractical except in the case of simply supported shells, the case that is derived herein. 
Therefore, an approximate theory, which takes into account general variations of displace- 
ments across the thickness, is developed for the static and dynamic analyses of laminated, 
composite, circular cylindrical shells with general boundary conditions. The approximate 
theory developed in this report is referred to as the "refined approximate theory." In the 
refined approximate theory the displacements are assumed to be piecewise linear across 
the thickness, and the three-dimensional problem is reduced to a two-dimensional problem 
in the circumferential and radial coordinates. Furthermore, a thin-shell theory is devel- 
oped without making any approximations in addition to the basic assumptions inherent in 
thin-shell theory. 

The accuracies of the refined approximate theory and the thin-shell theory are 
assessed by comparing results obtained by using them with the exact results for simply 
supported shells. 


SYMBOLS 

A A« A arbitrary constants used in particular part of refined approximate 

r’ 6 ’ z 

analysis 

a outer radius of shell, a*(p + 1) 

a*(i) inner radius of ith ply 
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a(i) 


a* (i) 
a 


B n (l,l), B 12 (l,l), . . . 

b 

b 


11 


C 66 


c 22 


c 66 


Dp, D n (l,l),D 12 (l,l), 


d r’ d 0’ d z 


F pr> Fqq, F q z , . . . 

fpr 5 f 60> *0z> • • • 

G(i,k) 


H(j,k) 


h*(i) 


stiffnesses of shell used in thin- shell analysis 
length of cylinder 

nodal distance or half -wavelength, b/n 
stiffness elastic constants 

modified elastic constants for thin- shell analysis 

stiffnesses of shell used in refined approximate analysis 

coefficients in power series for displacements in exact three- 
dimensional analysis 

integrals of stresses in refined approximate analysis 

integrals of stresses in thin- shell analysis 

arbitrary constants in expressions for displacements and 
stresses in exact three-dimensional analysis 

coefficient in power series expansion of displacements used in 
exact three-dimensional analysis 

thickness of ith ply 


hw = m 


h k, f, /3 indices 

M(/3) matrices occurring in particular part of refined approximate 

analysis 

m > n integers occurring in trigonometric expansions of displace- 

ments and stresses in circumferential and axial directions, 
respectively 
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ri7T3- 

b b 

0 null matrix 

p number of plies 

n number of layers into which middle ply of a three-ply laminate 

is split in refined approximate analysis 


Qee> • • • 


q rr (i), q 
q rz (D 


re 



q rr (2), q r0 (2), 

q rz (2) 


q 0r J q 00> q Qz 


q zr’ q z0’ q zz 

q rrf^’ q r0f^^’ q rzf^^’ 



r 


integrals of applied stresses on edges 

applied tractions on inner surface of shell in r, 9, and z 
directions, respectively 

applied tractions on outer surface of shell in r, 9, and 
z directions, respectively 

applied stresses on edges 9 = Constant in r, 0, and z 
directions, respectively 

applied stresses on edges z = Constant in r, 9, and z 
directions, respectively 

Fourier coefficients of applied pressures 


radial coordinate 


T 


kinetic energy 


t time 

U (i) difference between radial displacements at inner and outer 

r v 

surfaces of ith ply used in refined approximate analysis 

U^(i) difference between circumferential displacements at inner and 

outer surfaces of ith ply used in refined approximate analysis 
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U (i) difference between axial displacements at inner and outer 

z 

surfaces of ith ply used in refined approximate analysis 

Ur O ,U 0 o,U z o radial, circumferential, and axial displacements of inner 

surface in thin- shell analysis 

u r ,u^,u z radial, circumferential, and axial displacements, respectively 

V(l,k,i), V(2,k,i), . . . coefficients of differential operators in governing differential 

equations of refined approximate analysis 

W strain energy 

O 

W e work done by applied loading 

work done by applied tractions on inner and outer surfaces 
of shell 

work done by applied stresses on edges 
z axial coordinate 


(WJ 


e' surfaces 


^e^ edges 


Z = # 


a(k) 

r 0 (/3), r z (/3) 

n 0 (js), n z (/3) 

I 0 (/3), T z (j8) 

f z (/3) 


indicial constants in Frobenius method of solution of differential 
equations used in exact three-dimensional analysis 


matrices occurring in system of ordinary differential equations 
in complementary part of refined approximate analysis 
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Krone eke r delta 


6(j,i) 


e rr> e 99> e zz> 
e rz> e r 9’ e 9z 


strains 


£(j,k) =1 for k * j and £(j,k) = 


T? r (i,/3), V e (i,P), 
V z (h P) 


X 


R - a(j) 


h(j) 


for k = j 


functions of Z occurring in complementary part of refined 
approximate analysis 

circumferential angular coordinate 

nondimensional frequency parameter, equation (98) 




functions of 9 occurring in complementary part of refined 
approximate analysis 


p(i) 


mass density of ith layer 


a rr> 

a rz ? 


°99> 

a r9> 



t(/3) 

4 ) t } & 9’ *^z 


X r (i,k), X 0 (i,k), 
Xz(b k )> X rr (i,k), 

XflflCijk), X zz (i,k), ^ 
Xrz(b k )> Xr0(i?k), 

X0z(^ k ) 

Cl 


stresses 

coefficients occurring in summation for coefficients V(/3,k,i) 

functions of radial coordinate occurring in double trigonometric 
expansion of displacements used in exact three-dimensional 
analysis 


functions which are power series in radial coordinate occurring 
in double trigonometric expansion of displacements and 
stresses used in exact three-dimensional analysis 


frequency of vibration 
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A tilde (~) under a symbol denotes a matrix. Primes denote differentiated quan- 
tities. Numbers in parentheses are indices. 

STRESS- STRAIN RELATIONSHIP 

Most of the materials commonly encountered in structural practice are either iso- 
tropic or orthotropic. But in some cases during analysis there arises a need to choose a 
coordinate system which is different from the axes of orthotropy. In such cases the 
stress- strain relationships exhibit anisotropy. As shown in reference 6, the anisotropic 
stress- strain relationship that might be encountered in the analysis of fiber- reinforced 
components is of the type 



If the material is orthotropic, the elastic constants C2g, £35, an< ^ C45 are zero * 

STRAIN- DISPLACEMENT RELATIONS 
The strain- displacement relations in polar coordinates are (from ref. 7) 
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REFINED APPROXIMATE THEORY 


Development of Theory 


An approximate theory for the statics and dynamics of laminated, anisotropic, 
circular cylindrical shells (figs. 1 and 2) is developed. The displacements are assumed 
to be piecewise linear across the thickness; that is, the displacements of the jth ply 
are given by 


j 

u r (j) = a 21 U r (k) ?(j,k) 
k=0 


3 

U 0 G) = a 22 u eW f 

k=0 


u z (j) = a 


3 

22 u z (k> ?a,k) 


(3) 


where j = 1, 2, . . ., p and 


S(j,k) = < 


1 for k ^ 3 
R - a(j) 


h(j) 


for k = j 


In equations (3), U r (k), U^(k), and U z (k) are functions of 6 and z and are independent 

of r; U r (0), U 0 (O), and U z (0) are the displacements of the inner surface; U r (k), U 0 (k), 

and U (k) are the differences between the displacements of the inner and outer surfaces 
z 

of the kth ply at a given 6 and z. 

A variational approach is used to obtain the governing differential equations and 
boundary conditions. The variation in strain energy W s , the kinetic energy T, and 
the work done by the applied forces W 0 due to virtual displacements are calculated 
and substituted into the variational condition (known as Hamilton’s principle) 
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(4) 


rh 

(6T + 5W e - SW s )dt = 0 


where t Q and t 1 are the initial and final times, respectively. The strain energy W g 
is given by 


,a(j+l) 


W s = 


~ H f f f prrQ) e rr«) + W 3) + e zz (3) 

2 0=1 J e J z J a(j) 

+ CT rz (j) e rz (i) + CT re (3) e r0 (j)+ CT ez (j) e 0z (j)]RdRd0dZ 


From equations (2) and (3) the strains can be expressed as 


,.v U r® 
e rr(l) - 


3 

00®=! I] + U r (k) C (j,k) 

n L - 


;(j) = II 


8U z (k) , % 

— ^ ?(j,k) 


n U z® 

e rz (j) = + 

h(j) 


L ^ «),B 




Ufl(j) 


h L 


aU r (k) 


- u 0 (k) «j,k) 


3 r 

V au 0 (k) i au z (k) - (ik ^ 
l—> az r a e ' 

k=0 L J 
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Substituting for the strains from equations (6) and integrating with respect to R give 


W„ = 


a 


J J 

J 6 J z j=l k=0 


| F rr (j,3,l) U r (j) 6(j,k) + F re (j,3,l) U e (j) 6(j, k) 


+ F (j,3, 1) U (j) 6(j, k) + F eg 0,1,6^ 


^ * M" 


, „ aU 7 (k) , 4 

+ F ZZ G> 2 > 5 1> -^-+F re (j,l,5 1 ) 


au r (k) 

-^-u e (k) 


v au_(k) .. „ . , 8U„(k) 

+ F rz (j,2, 6i) — + Fg z (],2, Sj) aZ 




where 6^=1+ 6(j,k). The integrals F are integrals of stresses defined as 
(FgG.M); Fg(3,2,% 3,^)} 


(7) 



where 

j3 = rr, 99, . . ., 9z 

Stresses can be written in terms of displacements by use of equations (1) and (6). 

By substituting for stresses in equation (8) and integrating, the integrals F can be 
expressed in terms of displacements. (See table I.) 

Taking the variation of the strain energy W s , integrating by parts, and then changing 
the order of summation yield 


12 






The kinetic energy T is given by 


T = 


o P a(j+l) 

i j p<*> 

2 j-l J 6 z a(j) 


9u r (j)“ 

at 


+ 


" 9 u 0 (j)“ 

~~9t 


+ 


8U z (j)' 


•\~i 1 2 


at 


R dR ae dZ (10) 


Taking the variation gives 
6T = a 3 2] 


f a(j+l) 

J P(j) 

3=1 '0 z a(j) 


au r (j) a6u r (j) 
at at 


au 0 (j) a6u 0 (j) au (j) asu (j) 


+ 


at at 


at at 


R dR d0 dZ 


After integration by parts, 

t- t. 


r 


5T dt = -a" 


'0 


'0 


9 2 u 0 (j) 

at 2 


P 



3=1 

'*0 *^z 


2 


a u, 

(j) 

+ ; 

at 1 


i 


a(j+l) 


p(j) 


2 - « u z<» 


8 u r (j) 

^ - 5u (]) 

at^ 


R dR d0 dZ 


(ID 


Substituting equations (3) into equation (11) and integrating with respect to R give 


ti 


5T dt = -a" 


/* £ II 

k=o J e J z 


p 


E v 3 ’ 2 ’ 5 ^ 


U~k 


5U r (k) 



P 


P 

"\ 

+ 

E F pe«> 2 > 6 1> 

6U 0 (k) + 

E F PZ «> 2 - 6 i) 

5U z (k)^ 


J=k 


J=k 

J 


( 12 ) 


(See table I for F pr , F p0 , and F pz .) 
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The total work done by applied loads W e is the sum of the work done by applied 
tractions on the inner and outer surfaces of the shell (W e ) sur f a ces and the wor one 
by applied stresses at the edges (W e ) e dg e s> that 

W - (W e ) sur f aces + edges 


( 13 ) 


Now if the shell is subjected to dynamic tractions q rr (l), and q rz^ a 

q (2), q rfl (2), and q rz (2) on the inner (R = a(l)) and outer surfaces (R = ) o e s e 

in the r, 0, and z directions, respectively, then the virtual work done by these applied 
tractions is 


a(wj 


e y surfaces 


= a 



P 

Jj L 

J 0 J z /3-r,Q,z 

-q (l)a(l)8U„(0)+q r „(2) 2 Z 6 V k) 
P k=0 


d e dZ 


p 


= a 


f r 21 2Z [ q rp (2) ' q r(3 (1) 6(k ’ 0) a(1) ] 6U /3 (k) * dZ (14) 
J 9 Jz /3= r,0,z k=0 


If the edge z = Constant is subjected to stresses q zr , q Z 0> and q zz and the e g 
0 = Constant is subjected to stresses q 0r , q 00 , and q 0z in the r, 6 , and z irec 
tions, respectively, the virtual work done by these stresses is 


6(WJ 


e' edges 


P 

a(j+l) 

j 

V \ 

• s f L L 

J q ^ 

/ , 6U (3 (k) ?(i,k) 

J 9 j=l j3=r,0,z 

a(j) 

_k=0 


RdR d0 


P 


a(j+l) 

i 

\r — i 

** s r L 

E 

/ q ^ 

/ , 6U (3 (k) ?(j,k) 

4 j=i 

/3=r,0,z 

J a(j) 

k=0 


dR dZ 
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Integrating with respect to R and rearranging yield 


P P 




6 ^ W e^edges " a3 f Z_i “ u /3'‘ v ^z/3' 

J 0 ( 3 =r,e,z k=0 j=k 


H H 5 V k) q z/ 3 o.si) 


do 


P P 


+ a 


f H H i )dZ 

/3=r,0, z k=0 j=k 


where 


V j ’ 6 i } “ 


a(j+l) 


a(j) 


q z/3 R 


R - a(j) 

h(J) 


6 r l 


dR 


a(j+l) 


J a(j) 


R - a(j) 
h(j) 


V 1 


dR 


( 15 ) 


(16a) 


(16b) 


Equations (9), (12), (13), (14), and (15) are substituted into variational condition (4). Since 
the virtual displacements are arbitrary, the coefficients of 6U r (k), 6U$(k), and 5U z (k) 

(k = 0, 1, 2, . . p) can be equated to zero to yield the governing differential equations 
(3p + 3 in number) and boundary conditions (3p + 3 on a boundary). 

The governing differential equations are 
P 

H f- a2F pr (3 ’ 2 ’ 6 1 > - 6 ti> k) Fri-ti ’ 3 ’ 13 

j-k 

^66 5 1^ + ~S6 R r0^’* ,<5 l^ + ~dZ F rz^’ 6 1^ 

+ q rr ( 2 ) " q rr (l) 5 ( k >°) a(l) = 0 (17) 
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X r-a 2 F p0 (j,2,5 x ) - 6(3, k) F r0 (j, 3, 1) + F r0 (j, 1, 6{) 
j=k L 

+ '§$ F ee^’ 1,5 i > + ~§Z F ez^ , 2 ’ 5 i^J 


+ q r0 (2) - q re (l) 6(k,0) a(l) = 0 


( 18 ) 


X f-a 2 F pz (i, 2 > 6 l) ' 6(i ’ k) F rz (3 ’ 3>1) + M W 3 ’ 1 ’ 6 !* 


j=k 


+ ^ Fzzti - 2 ’ 6 ! 3 


+ q rz (2) - q rz (D 6<M) a W = 0 


for k = 0, 1, 2, . . • > P* 

The boundary conditions are as follows: 

For 6 — Constant, 


P * 

T. E Q 0r (3 > 5 l ) = 

j=k i =k 


0 or 6U r (k) - 0 


(19) 


(20 a) 


X FeeO.i.Oi)- E ( W 3 ’ 5 i ) = 0 or 6Ue(k) = ° 

j=k 3= k 


(20b) 


T F fl z (j,M x )- EQez«> 5 i )=0 0r 5 U z«=° 

j=k i= k 


(20 c) 
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For z = Constant, 


P 


22 F rz (j,2 )6l )- 21 Q zr a,5 1 ) = 0 or 5U r (k) = 0 

j=k j=k 


(21a) 


22 F 0 Z G’ 2 ’ 6 1>- 22 5 1> = 0 or 5Ue(k)=0 

j=k j=k 


(21b) 


j=k 


P 


L F z Z (5* 2 ’ 6 l)- H Q zz (j.«l) = 0 °r 6U z (k)=0 


(21c) 


j=k 


where k = 0, 1, 2, . . p. It is noted that for cylindrical panels, the boundary conditions 
for both Q = Constant and z = Constant are satisfied, whereas for cylinders, only 
the boundary conditions for z = Constant are satisfied. 

By substituting for the integrals F in terms of displacements, the governing differ- 
ential equations (17) to (19) can be written as 


P 


22 

i=0 


" V (1 , k, i) + V(2 , k, i) ^1+V(3,k,i) « 2 

dt* 


d6‘ 


3Z- 


+ V(4,k,i) -JL-g + V(5, k, i) 


U r (i) 


V(6, k, i) ~ + V(7,k,i) JL 


Ufl(i) 


+ 


V(8, k, i) -J— + V(9, k, i) JL 


u z (i)f = ^(k> °) a(l) - q rr (2) (22) 
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£ 

i=0 


V(10,k, i) -jjs + V(U,k,i) gZ 


U r (i) 


+ 


- V(1 , k, i) 4, + V(12,k,i)^ + V(X3,k,i) 


+ 


V M4 k i) a2 - + V(15,k,i) 
V(14,k,i; dQ dZ 


J 


U fl (i) 


V(16 


k i) _a^_ +v(17,k,i) + V(18,k,i) dQ qZ 

’ ’ d 6 2 92 


+ V(19,k,i) 


_J 


u (i) > = q r0 (l) 6( k >°) a W " 


(23) 


P f r a ' 

J V(20,k,i) +V(21,k,i) 

i=0 


U r (i) 


+ 


+ 


V(16 


k i) -§L + v(17,k,i) -2^ + V( 18 ,k,i) g 8 gZ +V(22,k,i) 

ae^ 9/1 


Un(i) 


-V(l,k,i) + V(23,k, l) ® 2 +V(24,k,i) 

at 


u 


+ 


V(25, k, i) + V ^ 26 ’ k> ^ 


U (i)f = 5 ^ k?0 ^ " q rz^ 


(24) 
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In equations (22) to (24) the V coefficients are pure functions of material properties and 
the inner and outer radii of the plies. (See table II.) They can also take into account var- 
iation of material properties within a ply. 

Procedure for Analysis 

The analysis of a cylinder or a cylindrical panel involves finding a solution of the 
governing differential equations (22) to (24) which satisfies the appropriate boundary con- 
ditions (eqs. (20) and (21)). The complete solution of the governing equations can be split 
into two parts: (1) a particular solution which takes care of the applied loading on the 
inner and outer surfaces and (2) a complementary solution which has arbitrary constants 
necessary for satisfying the edge boundary conditions. The procedure for obtaining the 
particular and complementary solutions is explained herein. 

Particular solution . - In some cases, depending on the applied loading on the inner 
and outer surfaces of the cylinder, the particular solution can be obtained by inspection 
of the governing differential equations. However, the general procedure is as follows: 

The applied loadings (amplitudes of loading in the case for forced vibrations) are 
expanded in a double Fourier series of the type 


OO OO 



. n7TZ 


n77Z 

b~ 


+ q rrf (/3,4) sin me cos 


n7rz 


b 


(25a) 


oO OO 


q r gO) = 2Z [" q r0f(& ^ sin me sln —— 

m=0 n=0 L b 


n7TZ 


+ q r0 f(£, 2) sin me cos . + q r0 fO, 3) cos me sin ^ 


+ q r 0f(& 4) cos me cos 


n?rz 


(25b) 
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rz 


(0) 


OO 00 

F JZ [ q rzf ((3 ’ 1) C0S ” 0 008 

m-0 n=0 


n7TZ 


+ q rzf O,2) cos m« Bui — 


+ q rzf (M) sin m ® C ° S 


n7TZ 

“IT 


+ q rzf (M) sinm0 sini ir] 


(25c) 


anri a t are functions of 

, 1 9 The Fourier load coefficients q rrf , <lr0f> an lrzI 

where /3 = 1, 2. loading some of them might be zero. 

m and n. Depending on the lo * ^ ^ also chosen f„ doubie Fourier senes 

displacements U r (i), z 


The 


as 


U r (i) = 


00 00 


L 1 

m=0 n=0 


F [A r (i, 1) cos m6 sln ^ 


„ _ nuz . a (i 3 ) sin m 9 sin n1TZ 

+ A (i, 2) cos m0 cos 
r ’ b 
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+ A r (i, 4) 


sin me cos 


— 

n77Z 

b _ 


(26a) 


o° 00 

U e (i) = £ E [A e (i,l)sinm 0 Sin 

m=0 n=0 L 


mrz 
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A (i, 2) sin m0 cos 
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n7TZ 
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A 0 (i,3) cos m0 sin 


n7rz 
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+ Ag(i) 4) 


cos me cos 



(26b) 
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u z (i) = 


OO oo 


L L 

m=0 n=0 


A (i, 1) cos m(9 cos 
z 


n7rz 

b 


+ A_(i, 2) cos m0 sin 
z 


n7rz 

b 


+ A (i, 3) sin m0 cos 
z 


n7rz 

___ 


+ A_(i, 4) sin m0 sin 
z 


mrz 

b 


(26c) 


where A r , A q, and A z are arbitrary constants. Equations (25) and (26) are substituted 
into the governing differential equations (22) to (24), and the various terms are grouped 

appropriately. Now, if the coefficients of cos m0 sin *! - — cos m.0 cos n ? z , 

sin m0 sin n7rz t and sin m0 cos are equated to zero in each of the equations, a 

b b 

set of 4(3p + 3) simultaneous linear algebraic equations is obtained for each m and n 
combination. The equations can be written in the matrix form 


M(l) g x 0 X J(4) 


A(l) 


K(l) 

0 1 M(2) g(3) 0 X 


A(2) 


K(2) 

0 1 J(2) M(3) 0 1 


A(3) 


5(3) 

J(l) 0 1 0 1 M(4) 


A(4) 


K(4) 


where 

AO) = {A r (0O), A r (l,0), . . ,A r (p,U); 

A 0 (O,/3), A 0 (l, /3), . . A 0 (p, /3); 

A z (0,P),A z (l,®, . . ., A z (p,/3)} 

K(/3) = ^q rr £(l> ^) a(l) - q rr f(2,/3)J, - <l rr f(2, /3) , . . 

|<lr0f(^> fi) a (^) - Qr0f^» > "^r0f^ J • • * > ^r0f^’ 

[q rzf (l,,3)a(l) - q rzf (2,|3)], -q rzf (2,/3), • • • > ' <l rz f ( 2 ’ 


(27) 


(28) 


(29) 


22 



Oi is a null square matrix of order 3p + 3, and M and J are square matrices of 
order 3p + 3. (See table III.) 

If all the plies are orthotropic, the J matrix becomes null, and correspondingly, 
equation (27) degenerates into four lower order equations, 

M(/3) A(/3) = K(/3) (30) 

where 0 = 1, 2, 3, and 4. 

Solution of the simultaneous algebraic equations (27) or (30) yields the values of 
A for each m and n combination. Summation of displacements and stresses with 
respect to m and n to desired accuracy completes the particular solution. 

Complementary solution. - For obtaining the complete complementary solution, 
the displacements U r (i), U$(i), and U z (i) are assumed to be in the form 


U (i) = 
r v ' 


GO 

J] r« r (i, 1) Sin 22^ + | r (i,2) cos M* 

n=0 L b 


+ sin md + ? 7 r (i, 2) cos mej 

m=0 


u e (i)= 2^ le 0 (i, 1) sin “2^. + i 0 (i, 2) cos 
n=0 *- 

OO 

+ 2] [v 1 ’ 1 * cos m 6 + 7Iq( i,2) sin mflj 

m=0 


U z (i) 


QO 

+ ^ [VzP’U Sin m0 + ^ z (i ’ 2) COS m0 ] 

m=0 


^ | z (i, 1) cos - n _— + £ z (i, 2) sin 

r\ b 


(31) 
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The £ are functions of 9 whereas the r\ are functions of Z. Equations (31) are 
substituted into the differential equations (22) to (24), now with a zero right-hand side, 

and the various terms are properly grouped. If the coefficients of sin n7rz cos n77Z , 

b b 

sin m0, and cos m0 are equated to zero in each of the equations, two sets of homogeneous 
simultaneous ordinary differential equations (6p +6 in number) are obtained. One set 
has 6 as the independent variable and | as the dependent variable. The second set 
has Z as the independent variable and rj as the dependent variable. The equations can 
be written in the following form: 


* 0 ( 1 ) n 0 (2) 

d 2 

5(D 

n 0 (i) te( 2 ) 

de 2 

5(2) 


r 0 (D i e (2) 

d 

«(D 

J 0 (D r 0 ( 2 ) 

de 

5(2) 


and 


t e (D 50(2) 



H 0 (l) f 0 (2) 


5(2) 


S Z (1) n z ( 2 ) 

d 2 

j(i) 

n z (i) t z (2) 

dZ 2 

2(2) 


r z (i) t z ( 2 ) 


vW 
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X z (i) r z (2) 

dZ 

2(2) 


+ 


f 2 ( 1 ) H z (2) 


2(1) 

S z (l) J z (2) 


2(2) 


(32) 


( 33 ) 
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where 


?o) = {i r (o,^), « r a,/3). • • • . ^(p./ 3 ); 
5 9 (0,/3), 5 e (l,/3), • • i e (p,/5); 
l z (0O), • • • > 

and 

7,0) = {i7 r (0,©, 7) r (l,/3), .... 77 r (P> ^ 
T/^Oj/ 3 )* 7,^(1, /3)'J * * * 7 IJgG?,/ 3 )* 
7J Z (0,^), 7, z (l,/3)> • • • . ^(P*/ 3 )} 


(34) 


(35) 


The matrices r 0 , T 0 , *0’ z ’ thp matrices 

matrices of order 3p + 3. (See table IV.) If all the plies are orthotropic the mat 

n e , Tg, Eg, n z , T z , and E z become null, and each of the equations (32) an ( 
reduces to two equations of order 3p + 3, that is, 


n 


T 


$ z , and H 


are square 


* e O) 5(/3) + r 0 O) ^ 50) + * g (0) m = 0 
de z ~ 


(36) 


^z(/3) 


^ tj(/3) + r z (/3) 

dZ^ ~ 


dZ 


v(P) + 


J z (/3) r?(/3) = 0 


(37) 


where ft = 1, 2. Solution of the homogeneous simultaneous ordinary differential equa- 
tions (32) and (33) or (36) and (37) yields expressions for 5 and t, 

One of the procedures of solving homogeneous simultaneous ordinary differential 
equations in closed form is given in appendix A. The solution of the differential equa- 
tions (32) and (33) can, in general, be written in the form 


|5 r (i>/ 3 )> Sgfbf 3 ); 


12p+12 

= 2Z {x r (i,ftl); x 0 (i,ftl); x z (i,ft*)} X U) 
i=i 


,g(j t)o 


(38) 
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and 


12p+12 

{vM,0)> v e d,P); v z (hP)} = 2Z 

X=1 


{y r (i,/3,X); y 0 (i,AX); y z (i,/3,X)} Y (X) e s ^ z 


(39) 


where X and Y are arbitrary constants. The functions x r ,x 0 ,x z and y r ,Y0,y z could 
be either constants or functions of 6 and z, respectively, depending on the multiplicity 
of roots g(f) and s(f), respectively. Equations (38) and (39) with equations (31) form 
the complementary solution. 

Further remarks on analysis . - In the previous sections the procedure for obtain- 
ing the particular and complementary parts of the solution to the governing differential 
equations was described for a general case. In many cases, depending on the loading, 
material properties, and boundary conditions, it might be sufficient to consider only a 
part of the general solution. After the appropriate part of the general solution is 
chosen, the boundary conditions are satisfied either exactly or approximately, as the 
case may be, to obtain the values of the arbitrary constants of the complementary part. 
Convergence studies can be carried out with respect to the number of terms chosen in 
the complementary part. 

Use of Refined Analysis with Simply Supported Cylinders 

In order to assess its accuracy, the present approximate theory is applied to the 
analysis of simply supported, laminated orthotropic cylinders (fig. 2), and the results are 
compared with those obtained with the exact three-dimensional analysis. 

The boundary conditions for a simply supported, laminated circular cylinder are 
on z = 0 and b 


^zzW = u r (i) = u #(i) = 0 (40) 

The boundary conditions (40) are identically satisfied by the first and third sets in 
the particular solution (eqs. (26)). If the loading on the cylinder is such that it can be 
expanded in the form of the first set in equations (25), then it is sufficient to consider only 
the first set in equations (26). Because the particular solution itself satisfies the support 
boundary conditions, the complementary part is unnecessary. Since the material is 
orthotropic, from equation (30) 

M(l) A(l) = K( 1) (41) 
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. .. , 0 _ 0 \ or forced-vibration (51 is known) problems, equation (41) is 

For static (51 - 0 or forced v t and displacem ents is summed 

solved to obtain the values of A. The series 


to the desired accuracy 
In free- vibration 


problems the right-hand side of equation (41) is zero; hence, 


(42) 


matrices, Af and A 2 51 2 , that 1S ’ 


(43) 


(44) 


M(l) A(l) = 0 

The matrix M(l) can be written as the sum of two 

M(l) = Af + Afjfl 2 

Substituting equation (43) into equation (42) gives 
Af A(l) + il 2 A 2 A(!) = ft 

Methods for numerical solution of characteristic matrix ® f each m and n 

are well known. Solution of equation (44) yields 3p + 3 frequencies 

combination. 

EXACT THREE-DIMENSIONAL ANALYSIS OF 
SIMPLY SUPPORTED SHELLS 

As mentioned previously, simply supported, 
cal shells (fig. — —rted 

analysis is possible. Th choosing the displacements and stresses in a 

. - . — rr- 

- “ rr.rirr * - ~ » — — 

conditions. 

Governing Differential Equations 

p she'll The three-dimensional equations of 
Consider a homogeneous orthotropic she! . 

equilibrium in a cylindrical polar coordinate system are (from ref. 7) 
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— — — ( r ^°r0) + — — <Jqq + — <Jn + p£2 2 u.0 = 0 

r 2 ar r& r dO 66 dz 6z P y 


(45b) 


i JL (ra rz ) + A A cr 0z + — c r zz + pft 2 u z = 0 

r ar r a0 az 


(45c) 


By use of stress- strain and strain-displacement relationships (1) and (2), the governing 
differential equations of three-dimensional elasticity in terms of displacements become 
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Solution of Differential Equations 


The displacements 

Up, 

uq, and u z are chosen to be in the form 
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where <fi r , <pQ, and 0 Z are pure functions of r. By substituting equations (49) in 
equations (46) to (48) and simplifying, a set of three homogeneous coupled ordinary differ- 
ential equations in which the independent variable is r is obtained; the equations are 
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The Frobenius method is now used to solve equations (50) to (52). The functions 0 r , 
(p 0 , and <p z are chosen in the form of a power series, that is, 


{< fir ; <p0> <#>z} = 22 R “ +j { H r^ H 0 (3) ; H z (j)} 

j=0 


(53) 


and are substituted into equations (50) to (52). This substitution results in 
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If the coefficients of R a 2 are equated to zero in equations (54) to (56), 
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For a nontrivial solution, the determinant of the coefficient matrix in equation (57) must 
be zero; as a result, the indicial equation 


(a 4 - 2VjO! 2 + v Q ) (C 44 a 2 - Cggin 2 ) = 0 


(58) 


is obtained where 


v 0 = 


C 22 (m 2 - l) 2 


11 


and 


v i = 


C 55^ C 22 + C ll^ + ™ 2 ( C ll C 22 ' 2C 12 C 55 ~ C 12^ 

2C 11 C 55 


Solution of equation (58) yields six roots for or, 


a(l) = ^v 1 + /vf- 


o;(2) = -o?(l) 


a(3) = ^ - -y/vf- 

o((4) = -Q?(3) 




o?(5) - m -y/Cgg/C^ 


ce(6) = -Q!(5) 


(59) 


Solving equation (57) for the constants H r (0), H^(0), and H z (0), one obtains (index k 

is added to denote the root) the following: 
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For roots, 1, 2, 3, and 4, 


~n 

H r (0,k) = G(k) 


Hg (0,k) = 


Cna(k) - Cggm - C 22 
m(Cgg + C 22 ) “ Iu(Cgg + C^2)“0 S ) 


G(k) 




H z (0,k) =0 


(60) 


and for roots 5 and 6, 


H r (0,k) =H 0 (O,k) =0 
H z (0,k) =G(k) 


(61) 


In equations (60) and (61), G(k) for k = 1, 2, . . .,6 are arbitrary constants. By con- 
tinuing the process of equating the coefficients of each power of R to zero (that is, 

R<*(k)-1 ( RO'(k), R a (k)+1 ; etc.) and solving the resulting equations, the constants 

H r (j,k), H 0 (j,k), and H z (j,k), for j = 1, 2, . . . and k = 1, 2, . . ., 6, can be 

expressed in terms of the arbitrary constant G(k) 

{H r (j,k); H e (j,k); H z (j,k)} = G(k) {d r (j,k); d 0 (j,k); d z (j,k)} (62) 

where d r , dg, and d z are functions of Q?(k), j, m, N, £2, and material properties. 
The functions d r , d$, and d z are obtained through recurrence relations (appendix B). 
From equations (53) and (62), 

6 °o 

b r ; <t>e-> <Pz] = 22 G(k) 2Z R “ (k)+i { d r G ’ k); d e (j ’ k); d z (j ’ k) } (63) 

k=l j=0 

If o?(k) happens to be a multiple root or an integer, the form of equation (63) changes 
slightly. The complete results are summarized in appendix B. 

If equations (49) and (63) are combined, the displacements u r , u 0 , and u z can 
be formally written (index i is added to denote the ith ply of the laminate) as 
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( 64 ) 


Expressions for x r , Xq, and x z are given in appendix B. By making use of stress 
strain and strain- displacement relations, stresses are obtained in the form 
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( 65 ) 


Expressions for X rr , Xqq, . . . are given in appendix B. 
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Boundary Conditions 


The displacements and stresses (eqs. (64) and (65)) satisfy the boundary conditions 


(eq. (40)) identically. The conditions that are to be satisfied are 
For the inner surface R = a(l), 

o'rr (!) = q rr (D, Orzfl) = q rz (D (66) 

For the outer surface R = 1, 

cr rr (p) = q rr (2), g t q (p) = q r0 (2), cx rz (p) = q rz (2) (67) 

For the interfaces R = a(f) where ft = 2, 3, . . p, 

o- rr (4 - 1) - cr rr (lO = - 1) - cr r0 (f) = cr rz U - 1) - cr rz (f) = 0 (68a) 

u r (£ - 1) - u r (f) = u 0 (£ - 1) - u 0 U) = u z (£ - 1) - u z U) = 0 (68b) 


Let the applied loadings be expandable in double Fourier series as in the first set of equa- 
tion (25). When conditions (66) to (68) are satisfied, a set of 6p simultaneous algebraic 
equations is obtained for each m and n combination. The equations can be written in 
the matrix form as 

PG = y (69) 

where 

G = {G(l,l), G(l,2) G(l,6); 

G(2,l), G(2,2), . . G(2,6); 


G(p, 1), G(p, 2), . . G(p,6)} (70) 

r = { q rrf( 1,:L ^ qref^’ 1 )’ qrzf^’ 1 ^ 23’ 

q rrf (2,D, q r0f (2,l), q rzf (2,l)} (71) 

(0g is a null column matrix of order 6p - 6) and 
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£(a(l),l) 
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°i ... 
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( 72 ) 


?2 


~2 


• L(a(p-l),p-2) -L(a(p-l),p-l) 0 2 

0 9 L{a(p),p-1) -L(a(p),p) 


2i 


m(i,p) 


In equation (72), 




X rr (i> 1)’ • • •> ^rr(^) 

Xp0(ij6) 

X rz (i, 1), X rz (i,2), . . Xrz^ 6 ) 


R=i^ 


(73) 


m(M) 


L(^i) 


X r (i, 1), X r (i,2), . . X r (i,6) 

X^(i, 1), X^(i>2), . . X^(i-j^) 

X z (i? !)> X z (i,2), . . X z (i,6) 


and 0-j^ and Og are null matrices of order (3 x 6) and (6 x 6). 


(74) 
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In static (£2=0) and forced-vibration (£2 is known) problems, equation (69) can be 
solved to obtain the constants G, and evaluation of equations (64) and (65) to the desired 
accuracy yields the displacements and stresses. 

In a free-vibration or wave- propagation problem, the right-hand side of equation (69) 
is zero. For a nontrivial solution, 

det P = 0 (75) 

The solution of this characteristic equation yields the frequencies of free vibration. 

For given circumferential and axial wave numbers m and n, there is an infinite 
spectrum of natural frequencies, each corresponding to a different thickness or radial 
mode. 

When either m or n equals zero, the characteristic determinant (of order 6p) 
degenerates into a product of two determinants of orders 4p and 2p. 

For m = 0, the determinant of order 4p corresponds to axisymmetric motion 

d the determinant of order 2p corresponds to torsional vibration 

. Furthermore, in the axisymmetric motion, if the nodal distance 

(b = b/n) is infinite, the radial and axial displacements become uncoupled, and the 
determinant of order 4p degenerates into a product of two determinants, each of 
order 2p. One corresponds to radial motion and the other to axial motion. 

For n = 0, the determinant of order 4p corresponds to plane strain motion 

, and the determinant of order 2p corresponds to thickness shear in 

the axial direction - = 0, u = un = 0 V 

W r / 

Procedure for Numerical Evaluation 

Static and forced-vibration problems . - If the geometric and material properties and 
the loading are known, the steps to evaluate stresses and displacements are as follows: 

(1) Express the loading in double Fourier series and obtain the Fourier load 
coefficients. 

(2) Compute the necessary x functions for radial coordinates R= a(l), . . ., a(p), 1. 

(3) Solve equation (69) to obtain the constants G(i, k), for k = 1, 2, . . .,6 and 
i = 1, 2, . . ., p. 

(4) Evaluate equations (64) and (65) to desired accuracy by summing the series with 
respect to m and n to obtain the displacements and stresses at points of 
interest. 
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Free vibration . - If m, n, and the geometric and material properties are known, 
the steps to obtain the frequencies of vibration to are as follows: 

(1) Assume a starting value for the frequency. 

(2) Compute the necessary X functions for R=a(l), . . a(p), 1. 

(3) Calculate the determinant | P | . 

(4) Adjust the frequency by a suitable amount and repeat steps (2) and (3) until the 
sign of the determinant for two successive values of frequency is different. 

(5) Linearly interpolate the frequency for which the determinant is zero, and by using 
this as the initial value and employing the well-known regula falsi technique, 
refine the approximate frequency to the desired accuracy. 

Some remarks on computation . - In all the problems — static, forced vibration, or 
free vibration - it is necessary to evaluate the X functions at the required radii, not 
only for solving equation (69), but also for calculating stresses and displacements at the 
required locations. In free vibrations the X functions will have to be evaluated for each 
trial frequency. The X functions, which are power series in the radial coordinate R, 
are slowly convergent. For the numerical results presented in this paper, the number of 
terms summed for good accuracy ranged between 80 and 140 depending on thickness, 
thinner cylinders requiring more terms. The numerical evaluation was time consuming. 

In free- vibration problems, in order to save computer time, frequencies were evaluated 
first by using the refined approximate theory, and these frequencies were then used as the 
starting values for solving the characteristic equation (eq. (75)). In contrast to the exact 
three-dimensional analysis, frequency calculations using the refined approximate theory 
and thin -shell theory took very little time. 


THIN-SHELL THEORY 


In this section the Fltlgge type analysis for homogeneous isotropic shells (ref. 8) is 
extended to anisotropic laminated circular cylindrical shells. The displacement distribu- 
tion consistent with the basic assumptions of thin- shell theory is 
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where U r Q, U^q, and U z q are the displacements at the inner surface of the shell. 
The strains are 
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The modified relationships between stresses cr zz , and a 0z and strains 

e 66> e zz> and € 0z f °r the 3 th ply are obt ained as follows: Equation (1) is solved for 

strains €pp* ^66* ^zz j and 
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Now, solving for stresses cr^, cr zz , and a 0z in terms of €qq, e zz , and e 0z and mak- 
ing use of the last three equations of equation (78) yield 
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A procedure identical to that for the refined approximate theory is now followed to 
derive the governing differential equations and boundary conditions. The strain energy 
W g is given by 
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Substituting for the strains from equations (77) yields 
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The stress integrals f can be written in terms of displacements. (See table V.) Taking 
the virtual variation of strain energy and then integrating by parts give 
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The virtual variation of kinetic energy, including the rotary inertia terms, is 
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where fp r , fp$, and fp Z are given in table V. The virtual work due to applied tractions 
on the inner and outer surfaces is 
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The virtual work due to applied stresses on the edges is 
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Equations (82) to (85) and (13) are substituted into equation (4), and the various terms are 
properly grouped. Since the virtual displacements are arbitrary, their coefficients are 
equated to zero to obtain three differential equations and four boundary conditions. The 
three differential equations are (after substituting for the integrals f in terms of 
displacement) 
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The quantities &22’ ®23’ * * • are defined in table V. The boundary conditions are 

as follows: 

For 0 = Constant, 


r 


di ee^ 

86 


az 


" f pr (2) + 


1 - a(l) 
a(l) _ 


= Q, 


aQ 


6 z 


6r 


az 


or 5U r0 - 0 





(90a) 

- Qqq 

or 

■C2*)- 

(90b) 

ffl 0 (l) = Q QQ 

or 

6 Ugo = 0 

(90c) 

%z( 2 ) = Q$z 

or 

5U z0 = 0 

(90d) 


For z - Constant, 


af 0z (3) 

86 


9 W 2 > . f 


az 


pr 


(3) + q rz (2) [1 - a(l)] 


= Q + or 5U n = 0 

zr 9 q rO 


(91a) 


W 2 > = Qzz °r 6 = 0 

W 1) = Q ze or 5U eo = ° 
f zz (D = Q zz or 6 U z0 = 0 
and at corners, 

% z (3) = Q ze +Q dz or 6 U r0 = 0 


(91b) 

(91c) 
(9 Id) 


(92) 


It is noted here that the thin -shell theory has only three differential equations in con- 
trast to the 3p + 3 equations of the refined approximate theory. Also, the number of 
boundary conditions that can be specified at any location on the boundary is 4 in thin-shell 
theory compared with 3p + 3 in the refined approximate theory. 
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Procedure for Analysis 


The analysis of a circular cylindrical shell now involves finding a solution of the 
governing differential equations (eqs. (87) to (89)) which satisfies the appropriate boundary 
conditions (eqs. (90) and (91)). The complete solution of the governing differential equa- 
tions can be split into a particular part and a complementary part. The procedures for 
obtaining the particular and complementary parts are similar to those for the refined 
approximate analysis. 

Use of Thin- Shell Theory With Simply Supported Cylinders 

In order to assess the accuracy of the thin- shell theory, it is applied to the analysis 
of harmonic vibrations of simply supported, laminated orthotropic cylinders (fig. 2), and 
the results are compared with those obtained by using an exact three-dimensional analysis. 
The displacements U r Q> Ugo> anc * U z o are assumed to be of the form 


II 

o 

u 

oo 

L 

OO 

L 


m=0 

P 

II 

O 

u e o = 

OO 

L 

OO 

L 


m=0 

n=0 

u z o = 

OO 

L 

m=0 

CO 

L 

n=0 


S cos m0 sin 
r 


n7rz 

b 




sin m0 sin 


n7TZ 

b 


> 


cos m0 cos 


n7rz 

b 




(93) 


If the applied loading can be expressed in the form of the first set in equation (25), then 
by substituting equation (93) along with the first set in equation (25) into equations (87) 
to (89), a set of three simultaneous algebraic equations is obtained for each combination 
of m and n. These equations can be written in the form 


Pl S +S2 2 P 2 S = CO 


(94) 


where 



(95) 
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S? = |( a(1) q rrfW - q rrf (2) - m q ref (2) + [1 - a(l)] Nq rzf (2)j ; 

( a(1) q rM (1) - ( E(1) q rzf (1) ' 9 r Z f( 2 >)j ( 96 ) 

and and P 2 are square matrices of order 3. 

The matrices Pj and P 2 are of the form 



where 

K 1 = -B 22 (M)m 4 - 2B 22 (2,2)m 2 - [2B 2 g(5,3) + Bg 6 (12,3)] m 2 N 2 - B 22 (3, 1) 
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k 6 = -B330, 1) N 2 - B 66 (3,l)m 2 

k 7 = Bp (9, 1) + Bp (7, 3) m 2 + B p (9, 3) N 2 


K 8 = Bp (10, 2) m 
Kg — -Bp (9, 2) N 

K 10 ~ 


k 1 1 = Bp (9, 1) 


See table V for B relations. 


in the case of static or forced-vibration problems, equation (94) * 


c q o nf i S The series for stresses and displacements are 

obtain the constants S , S 0 , and b z . me series 


summed to the desired accuracy. 

In the case of free vibrations, equation (94) becomes 


P x S + ft 2 ?2§ “ ft 


(97) 


Phis characteristic matrix equation is solved to obtain the frequencies of free Ration. 

and P 2 matrices is 3, the number of frequencies is also 3 for 


Since the order of 
a given m and n combination. 


NUMERICAL RESULTS AND DISCUSSION 


In this section the numerical results for the free vibrations of simply supported, 
laminated orthotropic cylinders are presented for two typical three-ply cylinders. 

R Tsult from the present refined approximate theory and the thin-shell theory are compared 
wUh results from the exact three-dimensional analysis. Since the simple- support condi- 
tions simulate the conditions at nodes in wave propagation in infinite cylinders, ' P re “ n 
results are applicable to the problem of wave propagation. In wave-propagation problems, 
the nodal distance b is equal to one-half the wavelength, and the wave velocity is equal 

to b Q/tt. 
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The frequencies are presented in terms of a dimensionless frequency parameter x, 
defined as 
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The first three frequencies of free vibration are presented in figure 3 and table VI 
for a thin cylinder (total thickness is 5 percent of outer radius) and a thick cylinder (total 
thickness is 20 percent of outer radius). The material and geometric properties are 
given in table VI. In both cylinders the middle ply is thicker and also of lower elastic 
moduli than the other two plies. In figure 3 the frequency parameter is plotted against 
the circumferential wave number m for various ratios of nodal distance to outer radius. 
In figures 4 and 5 the displacement distributions across the thickness are plotted for the 
thick cylinder. 

On attempting to classify the first three radial or thickness modes, it was found that 
in most cases not one of the three displacements was clearly dominating. In this respect 
the laminated composite cylinders differ from homogeneous isotropic cylinders, where 
the three thickness modes can be distinguished reasonably well as modes associated with 
either large radial displacements, large axial displacements, or large circumferential 
displacements. In the numerical results presented, for m = 0, the first mode was 
axisymmetric, the second mode was torsional, and the third mode was axisymmetric. 

In, contrast, when m = 0 in homogeneous isotropic cylinders, the first mode is generally 
torsional. 

As in homogeneous isotropic cylinders, as the nodal distance approaches infinity, the 
first and second frequencies vanish when m — 0, and only the first frequency vanishes 
when m = 1. For m ^ 2, all the frequencies approach nonzero finite values as the nodal 
distance approaches infinity. The frequencies are not always monotonic functions of circun 
ferential wave number m; the nature of variation of frequencies with m depends on nodi 
distance and material properties. (See fig. 3 and table VI.) 

In axisymmetric modes (fig. 4), even when the cylinder is thick and the nodal distance 
is short, the true u r distribution across the thickness is nearly piecewise linear in the 
first two modes. Also, the true u z distribution is piecewise linear in the first mode 
and becomes nonlinear in the second mode. In the torsional mode the true u^ distributio] 
is nonlinear. 


50 



Exact 

Thin- shell theory 



o 1 


2 3 4 



2 3 4 


2 3 4 


m „ Third frequency 

First frequency Second frequency 

(a) Thin cylinder: Total thickness, 5 percent of outer radius; a(l) = 0.95; 

a(2) = 0.955; a(3) = 0.995. 

j. ,, function of circumferential wave number m 

Figure 3.- Frequency parameter X as a function o 

for various nodal distances. (See table VI for material properties.) 


Fii?s"b fz*GQUGncy SGcond. fi^Gc^iiGncy Thix*cL fT'GcjiiGricy 

(b) Thick cylinder: Total thickness, 20 percent of outer radius; a(l) = 0.8; 

a(2) = 0.82; a(3) = 0.98. 


Figure 3.- Concluded. 





Exact 

Thin shell 


Refined approximate theory 

Refined approximate theory with 

mid- layer split into two 
layers 

Refined approximate theory with 

mid-layer split into four 
layers 



Second mode 



Figure 4.- Displacement distribution across thickness in a thick three-ply cylinder, 
m = 0; b/a = 1. (First and third modes are axisymmetric, u 0 = 0. Second mode 
is torsional (u r = u z = 0). Displacements have been normalized with respect to their 
maximum values.) 


When m = 1 (fig. 5), the true displacement distributions across the thickness are 
nearly piecewise linear for the first mode even when the cylinder is thick and the nodal 
distance is short. But for the second and third modes, especially when the cylinder is 
thick, the true displacement distribution across the thickness within each ply is nonlinear. 

The thin-shell theory overestimates frequencies (fig. 3 and table VI). Table VII 
summarizes the influence of thickness, wave number m, and nodal distance on the 
accuracy of frequencies obtained from thin-shell theory. In general, the error in the 
second frequency is less than the error in the first and third frequencies. The thickness- 
wise displacement distributions predicted by thin-shell theory are highly erroneous 
(figs. 4 and 5). 

The first frequency obtained from the refined approximate theory is very close to the 
exact value, even for thick cylinders, for all values of m, nodal distance, and material 
properties (table VI). The second and third frequencies are also reasonably close to the 
exact values. For the thick cylinder, the maximum errors in the second and third fre- 
quencies obtained from the refined approximate theory are 13 percent and 6 percent, 
respectively. In comparison, the corresponding maximum errors in the second and third 
frequencies obtained by thin-shell theory are 60 percent and 165 percent. The accuracy 
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Exact 

Refined approximate theory with 

mid-layer split into two 
layers 

Refined approximate theory 

Thin- shell theory 



(a) Radial displacement. 


Third mode 




(b) Circumferential displacement. 


Figure 5.- Displacement distribution across thickness in a thick 
three-ply cylinder, m = 4; b/a = 1. (See table VI for geo- 
metric and material properties. Displacements have been 
normalized with respect to their maximum values.) 
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Refined approximate theory with 

mid-layer split into two 
layers 

Refined approximate theory 

Thin- shell theory 



Figure 5.- Concluded. 


. . bv the refined approximate theory can be further improved by 

; the frequencies obtained by „ Ues . The second and third frequency 

putting the thicker ply into a number o q£ equal thickne ss are given in 

alues obtained by splitting the middle ply in ^ xima te theory can be brought 

.b.. ««■ M u,nei« lh " omb „ of plie, *«> “» 

,s close to exact values as desired by ^creasing J^’ccuracy with increase in P m is 
niddle ply is split (fig. 6). The improvemen { computa tion increases slightly 

■apid at first but slows down subsequently. The amount 

vith increase in P m * . , thpnrv for a given m 

The number of frequences predieted byjh, including the artificial 

refined approximate theory correspo the0 rv might identify some modes which 

The refined dLpl.eement di.iribution 

- 

thicker plies into thinner plies (figs. 4 and 5). 
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(a) Second frequency. (b) Third frequency. 

Figure 6.- Percent error in the frequency parameter as a function 
of the number of layers p m into which middle ply is split, 
b/a = 1; thick cylinder: Total thickness, 20 percent of outer 
radius; a(l) = 0.8; a(2) = 0.82; a(3) = 0.98. (See table VI 
for material properties.) 

CONCLUDING REMARKS 

A refined approximate theory for the static and dynamic analysis of finite laminated 
composite shells was developed. The analysis was reduced to a two-dimensional problem 
in the axial and circumferential coordinates by assuming piecewise linearity of displace- 
ments across the thickness. The governing differential equations and the boundary condi- 
tions were derived by using a variational approach. General solutions of the governing 
differential equations were developed in trigonometric series form. Analysis of finite 
laminated shells with general boundary conditions now involved satisfying the boundary 
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. , ^ OT , + p-eneral solution. This theory 

conditions by making use of the appropna ® Qn a bounda ry in contrast to only 

allows 3 + (3 x Number of plies) boun zry ^ ^ re fined approximate theory 

ion o me boundary 'conditions and, in turn, the results close to 
boundaries should be obtainable to a high degree of accl ™^ ^ ^ developed> 

An exact three-dimensional analysis of Simpl J * U ^ P ° a glng the accuracy of the 
and the results from it for free The thin-shell theory developed 

” f^rTZTr^ZZZ an^T approximations beyond the basic assumptions inherent 

short nodal distances and high circumferen (thickness about 5 percent of 

theory was found to be accurate for wave numbers, 

the radius) only in a narrow range of nodal distances an ^ ^ can be 

The accuracy of the results obtained from th ® re “ e .^ o & number oJ thinner plies, 
further improved, if needed, by splitting the jac « P ^ ^ ^ number of plies 

The accuracy can be studied throug *“ pr refined approximate theory can be 

S^=ems of laminated composite, circular 

SCI! sheTs. y The theory can be extended to arbitrary shells. 

Hampton, Va., August 20, 1973. 
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APPENDIX A 


SOLUTION OF HOMOGENEOUS SIMULTANEOUS 
ORDINARY DIFFERENTIAL EQUATIONS 


A procedure for solving a system of homogeneous simultaneous ordinary differential 
equations of the type of equations (32) and (33) is developed herein. The system of second- 
order differential equations is converted to a system of first-order equations, which can 
then be solved by following standard techniques. Consider a system of a equations, 


A — $ + bA$+C$ = 0 
~ dy 2 ~ ~ d y ~ 


(Al) 


where y is the independent variable, is the column matrix of dependent variables, 
and A, B, and C are the coefficient square matrices of order a. Now by choosing 


(p= — 3 > 
~ dy ~ 


(A2; 


equation (Al) can be written in the form 


aA. 0+B0+C$ = O 
- dy ~ 

Equation (A3) can be written after some manipulation as 


(A3 


-A -1 BcJ - A -1 C$ = (f) 

Combining equations (A2) and (A4) yields 


(A4 


where 



d 




_ $ 

dy 

X 



-A -1 C 


-A -1 B 


(A 5 


(At 
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APPENDIX A - Concluded 

a null square matrix of 

r \a\ T is a unit matrix of order ot, and 

" VJn ol m «» « 

(J; £} = {J; efXy 

„ / * n\ lnto equation (A5) and simplifying give 
Substituting equation (A7) into q 


(A7) 




= M- ^ 


D | 

L7_ 


L - J 


(A8) 


Solution of characteristic matrix equations such as 

„ 1. 2 «,«her.win» *> v1«,o. 

can be finally written in the form 


equation (A8) is well known. Since the 
|i. Thus, the solution of equation (A ) 


<| = 


£ S(f) T(t) e* {l)y 
1=1 


(A9) 


where E(£) 
constant. 

In cases 
roots can be 


is the eigenvector corresponding to root 


(iti), and T(t) is the arbitrary 


where multiple roots occur, the special solutions corr 
obtained by following the usual procedure. 


esponding to such 
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APPENDIX B 


X FUNCTIONS OF EQUATIONS (64) AND (65) 

In this appendix, expressions for the X functions occurring in equations (64) and (65) 
are given. The form of the solution given by equation (63) is true if the roots are nonzero, 
are nonmultiple, and do not differ by an integer. When any of these conditions are not 
satisfied, as in the cases discussed later, the form of solution is different from that of 
equation (63), and special procedures must be followed to obtain the solution. 


Ordinary Case 

When equations (63) and (64) are compared, it is seen that the displacement functions 
X r , Xq, and X z are 

OO 

|x r (k); X 0 (k); X z (k)} = E R °^ + ^ {d r (j, k); d 0 (j, k); d z (j,k)| (Bl) 

3=0 

The functions X rr (k), X 00 (k), . . . can be written in the form 


where 


Xrr(k) ~ E [^dj.&k) + i// 2 d 0 (j,k) + ^ 3 d z (j,k)] 

3=0 


(B2) 



^1 

^2 

^3 

Xrr(k) 

C u [a?(k) + j] + C 12 

C 12 m 

-C 13 NR 

x 00 ( k ) 

C 12 [«(k) +5] +C 22 

C 22 m 

-c 23 nr 

x zz W 

c 13 [“(k) + J] +c 23 

C 23 m 

-c 33 nr 

X 

*■* 

IS 

c 44 NR 

0 

C 44 [o?(k) + j] 

X r 0( k ) 

- C 55 m 

c 55 [«( k ) +3-1] 

0 

x 0z ( k ) 

0 

C 66 NR 

' C 66 m 


and d r , dg, and d z are given in table VIII. 
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APPENDIX B - Continued 

Special Case: m = 1 

onri the original solutions, given in 

When m = 1, roots «(3) and ^LeTwo indicial constants become identical, 

equations (Bl) and (B2), correspon 1 taine d but for the indicial constant oi , a 

For the root *3), the original solution is reU ^ for ^ is found . The 

special solution, which is linearly in e P e original solutions (eqs. (Bl) and (B2)) 

With respect to a, then substituting k equa tion, the index k is retained 

thus obtained is same equation will be used later with a differen 

instead of replacing it with 4 

value for k)z 
With k - 4, 


oO 

jXr(k); Xq{ k); X z (k)^- - ^ {]- r 


, k) log e R + 


d^(i,k)]; 


[d 0 (i,k) log e R +dg(i,k)]; [d z (S,k) log e R +d z«’ k) ]} 


(B3) 


Functions Xrj*(k), X 00 (k), 


can be written in the form 


X r rW = 


oO 

z 

j=0 


r )- 1 [^{d r (j,k) log e R + d^(j,k)} + log e R 


+ d 0 (i,k)j + ^ 3 { d z ( j> k) l0geR + A ' zih ^ 
+ 4/ 4 d r (i,k) +^ 5 d 0 (j,k) + ^ 6 d z (i ’ k) j 


(B4) 
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APPENDIX B - Continued 


where 



^4 

^5 

^6 

Xpi*(k) 

c n 

0 

0 

Xee(k) 

C 12 

0 

0 

Xzz(k) 

C 13 

0 

0 

Xrz(k) 

0 

0 

C 44 

Xrfi(k) 

0 

C 55 

0 

X$z(k) 

0 

0 

0 


In equations (B3) and (B4) the primes denote that the value has been differentiated with 
respect to a and then the limit has been taken as a - 0. The functions d r , d e , an 
are given in table VIII. 

Special Case: m = 0 

When m = 0, o<3) = 1, a(4) = -1, and o<5) = o<6) = 0. For m = 0, that is ^ = 0, 
displacements u r and u z are not coupled with u e ; also, the differential equations (50) 
and (52) are not coupled with equation (51). The motion that is associated with u r and 
u 7 is axisymmetric, and the associated indicial constants are o(l), a(2), o<5), and <*6) 

The motion associated with u 0 is torsional, and the corresponding indicial constants are 

a(3) and o?(4). 

Axisymme tric motion .- Among the four associated indicial constants, o<5) and o<6) 
are zero. Therefore, the original solutions, given in equations (Bl) and (B2), correspon - 
ing to these two indicial constants, become identical. For o<5) the original solution 
(eqs. (Bl) and (B2)) is retained, and a special solution corresponding to o<6) is found. 
This solution is obtained by differentiating the original solution with respect to a and 

then substituting k = 6 and o(k) - 0. The solution is the same as given in equa- 

tions (B3) and (B4), now with k = 6. 

Torsional motion.- The two indicial constants corresponding to torsional motion are 
o(3) and o(4), which are equal to +1 and -1. Since the indicial constants differ by an 
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APPENDIX B - Continued 

Integer, it would be ?L S "a ' 

the general solution, in this case it is -u Qn for torsio nal vibration is (from 

equation and solve it directly. The dill 

eq. (5D) 


'55 


dR 


^55 + (pSl2 a 2-C 66 N 2 ) +C 55 i ^ 

R 2 


<#>0 = 0 


(B5) 


Rearranging gives 


r ft 2 + R -i- + (R 2 V - 1) 

* 




$9 


= 0 


(B6) 


where 


R = ZTR 

uJ = -/Tw' 

V = 1, "1 




( 0)2 > 0 , a) 2 < 0 ) 




J 


(B7) 


and 


2 pH 2 a 2 _ C 66 n 2 

to =^~- 7 ; 

C 55 C 55 


(B8) 


2 > 0 and a modified Bessel equation when 
Equation (B.) is a Bessel equ^n 1 when . >0 ^ ^ glven by 

2 < 0 . The two solutions of equation , 


Xg(3) 


Ji(R) 


I,(R) 


(to 2 > 0) 
(to 2 < 0) 


(B9) 
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APPENDIX B - Continued 


X 0 ( 4 ) - 



(0,2 > o) 

(o, 2 < 0) 


(BIO) 


where and are Bessel functions of first and second kind of order 1, and 1^ and 

are modified Bessel functions of first and second kind of order 1. 

Among the other X functions, X r , X z , X rr , Xqq, X zz , aa d ^rz are zero, and 


x r e(3) - c 55 co . 

x r0 (4) = c 55 ut[^5) 
X ez (3) = CggNJ^R) 
X 6z ( 4) = CgeNY^R) 



(Bll) 


When o, 2 < 0 in equations (Bll), J is replaced by I and Y by K. 

Special Case: m -^Cgg/ = Integer 

When m -JC qq/C^ is an integer, the indicial constants o(5) and o(6) differ by an 
even integer. In this case it would be impossible to calculate the solution corresponding 
to of(6) by using equations (Bl) and (B2). The special solution is obtained by multiplying 
the original solution (eqs. (Bl) and (B2)) by a - a(6), differentiating with respect to a, 
and then substituting oi - o;(6). Now, 

i3-i 

{X r (6); X e (6); X z (6)} = R°( 6)+ j {d r (j, 6); d e (j, 6); d z (j, 6)} 

j=o 

oo 

+ ^ R^ 6 ^ +;j |[d r (j, 6) log e R + di.(j, 6)J ; 
j=/3 

[d 0 (j, 6) log e R + dgQ, 6)] ; [d z (j, 6) log e R + d z (j, 6)] j (B12) 
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APPENDIX B - Concluded 

r ’ , r.\ ta\ pan be written in the form. 

where (3 = X(5) - X(6) = 2m V C 66 /C 44' x rr( 6 )> x re(6), • • • 


( 3 - 1 -I 

X rr (6) = 2^ R^ 6 ^" 1 jVidj.0,6) + ^2 d 0^’ 6 ^ + ^3 d z^j 6 )J 

j^O 


+ 22 R 0 ^ 6 ^ -1 |^4d r (j,6) + i//5d 0 (i,6) + i^ed z (h 

3=13 


6 ) 


(B13) 


+ [d r (j, 6 ) lo Se R + ^r(3> 6 )] 

+ [a e (j, 6) log e R + 4(j, 6)] + ^ [4(3, 6) log e R + 4(5, 6)]} 

In equations (Bl2) and (B13), 

d r (j, 6) = lim {[a - o(6)] d r (],6)j 

OL-~oi 6) 

4(3, 6) = & {[« - <«>] Vi> 6 )} at “ = 0(6) 

and d r , di-, . . . are given in table VIII. 

Special Case: Isotropic Material 

equations (Bl) and (B2). •> pan be obtained very 

solution to ,h, governing dM.r.nM «,n,.ion. («) to 48» ~nbe 

easily by tallowing the displaeement potential techm,ue (see ret., and ). 
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TABLE I.- EQUATIONS FOE F 


] 

The F (j 3 1) F r( o(;j,3, 1), . . . can be expressed in the form F rr (;j,3,l) = ^ Q< r (i) U r (i) + Kg(i) Ug(i) + K z (i) U z (i)] 
rr ’ ’ i=0 


F rr (j,3,l) 

F r 0 (j, 3,1) 

F rz (j,3,l) 

F r0 (j,l,6i) 

F 0e (j,l,6l) 

F 0z (3,l,6l) 

F rz (j,2,6i) 

F0 Z (j >^>®l) 
F zz a,2,6i) 
F pr (j,2,6!) 
FpeQ^.fi!) 
F pz (j,2,6i) 


K r (i) 


6(3,1)044(3, 7,1) + Di 2 (i,5,6 4 ) 
055(3,5,64) — + 045(3,4,64) — 


045(3,5,64) — + 044(3,4,64) 


BZ 

8 

9Z 


055(3,3,65) ~ + D 45 a,l,6 5 ) 
6(j,i) D 12 (j,5,6i) + 022(3,3,65) 
6(3,1)016(3,5,64) + 025(3,3,65) 
D 45(i. 1 . 6 5 )“+ D 44(l.2.«5 )^ 
6(3, i) 045(3,4, 64) +D 2 6(3,1,65) 

6(3,1)043(3,4,64) + 023(3,1,65) 

g2 

Dp(j,2,6 5 ) -2 

p at 2 

0 

0 


K e(i) 

— — — 

K z (l) 

042(1,5,64) + 045(3,4,64) • 

016(3,5,64) Jg + Dl3(j ^,64) ~ 

6(3,1)055(3,7,1) -055(3,5,54) 

6(j,i) 045(3,7,1) 

6(3,1)045(3,7,1) - 045(3,5,64) 

6(3, i) D 44 (3,7,l) 

6(3,1)055(3,5,64) - 055(3,3,65) 

6(j, 1) 045(3,5,64) 

022(3,5,65) -QQ + 026(3,1,65) 

026(3,2,65) -gg + 023(j,l,65) £2 

0260,5,65) -gg + 055(3,1,65) Q2 

066(3,3,65) — + D36(j,l,65) ^ 

6(3,1)045(3,4,61) - D45(j,l,65) 

6(j,i) 044(3,4,64) 

O260, 1 , 65) + 056(3,2,65) 

D 66(3,M5) Jg + D 36(3,2,65) 

D 2 3(3,!,65) Jg + D 36(3, 2 ,65) 

035(3,1,65) jjj + 033(3,2,65) — 

0 

0 

8 2 

Dp(j,2, 65) 

H 8t 2 

0 

0 

D 0 (j ,2,65) 

p 8t 2 


Note: 64 = 1 + 6(3, i); 65 = 64 + 6(j,i) 

The stiffnesses D are pure functions of the material properties and the inner and outer radii of the individual plies; 
that is, 

D ? (j,7,«] - £° +1> fi; Ri b 3:- 3^; -^3 c «<i> 


b(j) 


J a (j) l ’ R 5 h(j) ; h(j) S Rh (j)’ h(j)2/ "P 

where /3 = 11, 12, . . ., 66, p. When j3 = p, = p. If the material properties are constant across a ply, 

5 /3 


0^(1, 2,£) = C (3 (j)h(j) 
D^(j, 3,1) = C0(j) log, 


"l + h(3)ll 

L soil i 


D fl (j,l,£) = C ^— 0^(3, 3,2) = Cpfo) jl - ^ log e 

T 1 ] 


Oft(i,l,£) 

=-W“ 


h(j) 
\J + 1 


1 + 


Mj) 

a(j) 


0^3(3, 4, H) 


Dgfj, 2, £) 


h(j) 
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TABLE II.- COEFFICIENTS V 
P 

V(j3,k,i) = £ A(j,i) t(j8) ; 
j=E 

k = 1 for k = 0 , ic = k for k > 0 ; 

A(j,i) =0 for i > j, A(j,i) =1 for i ^ j 


D p (j,2,5 5 ) a2 

D55(i> 3 » 6 5) 

044 ( 3 , 2 , 65 ) 

2045(3,1,65) 

-6(j,k) 5(j,i)DiiO,7,l) - 60,k)D 12 0,5,64) - 6(j,i) D 12 0,5,6l) - D 22 (j,3,65) 
-60, k) D4 2 (j,5,64) - D 22 (j,3,65) - 055(3,3,65) + 6(3,1)0550,5,64) 
-60,k)D 16 (j,4,64) - 045(3,1,65) + 6(3,1)045(3,4,61) - D 26 (j,l,6 5 ) 

-6(3, k) 0450,5,64) + 6(3, i) D 4 50,5, 61) - D 26 (j,3,65) 

-6(j,k) 0430,4,64) - 023(3,1,65) + 6(j,i)D 4 4(j,4,6i) 

-6(3, k) 0550,5,64) + 055(3,3,65) + 6(3, i) D4 2 0,5,6 i) + D 22 (j,3,65) 
0450,1,65) + 60,i)D 16 0,4,6i) - 6(3, k) 045(3,4,64) + D 26 (j,l,6 5 ) 


D22OA65) 


066 ( 1 - 2 , 65 ) 

2D 2 60,f ,65) 

6(3, k) 055(3,5,64) - 6(j,k)60,i)D 5 50,7,l) - 055(3,3,65) + 6(j,i) 0550,5,64) 

d 260> 3 >65) 

036 ( 02 , 65 ) 

023(01,65) + D 6 60,l,6 5 ) 

-6(j,k)60,i)D4 5 0,7,l) + 6(3,1)045(3,5,6!) 

O26(0 3 >65) - 6(3, k) 045(3,5,64) + 60,0Di6(OMi) 

6(3,1)0430,4,64) - 6(j,k) D 4 4(j,4,64) + D 23 0,l,6 5 ) 

-6(3, k) 6(j,0D 4 5(j,7,l) + 6(j,k) 045(3,5,64) 

D 66(3,3,65) 

033(3,2,65) 

2036(3,1,65) 


26 - 6 (j,k) 6(j,0 D 44 0,7,l) 

Notes: 64 = 1 + 6 (j,k), 65 = 1 + 6(j,0 + 6 (j,k), 64 = 1 + 60,0, and the stiffnesses D are 

given in table I. 



TABLE in.- MATRICES MO) AND JO) 

Matrices MO) and 3 ( 0 ) 

intheform r Si o, ia u» .,<*! I 

M(/3) = ^ 4 0) ~6 ( ^ 

j^O) S 8 (® 

^ /ox * Q (/3) are given below in Which 

tprms of matrices *i(0)> • • *’ W 9 W 


General terms of 
k refers to row 


and i to column. 


MO) — — — — — ~ ~~ 

— — — -V(4,k,i) mN6 2 (/3) 63 O) 

V(l,k,l) - V(2,k,l) m2 - V(3,k,i) N 2 + V(».M 


V( 6 ,k,i) mS 2 0) 

-V(9,k,i) N5 3 (/3) 

-V(10,k,i) m5 2 (/3) 

V (l,k,i)Sl2 . V(12,k,i)m2 . V(13,k,i)N 2 + V(15,k,i) 
V(18,k,i) mN5 2 0) S 3 O) 

V(21,k,i) N 63 O) 

V(18,k,i) mN8 2 0) 63 O) 

V(l,k,i)fi 2 - V(23,k,i)m 2 - V(2 4,k,i)N 2 + V(26,M ^ 


V(7,k,i)N6 3 0) 

-V(8,k,i) m6 2 0) 

V(ll,k,i) N5 3 0) 

V(14,k,i) mN6 2 0) 5 ^( 0 ) 

-V(16,k.l)m2 - V(n,k.i)N 2 + V(19.k,i) 

-V(20,k,i) m5 2 (/3) 

-V(16,k,i) m 2 - V(17,k,i)N 2 + V(22,k,i) 
V(25,k,i) mN6 2 0) S 3 O) 


Not'e: H W - ^ 2 ( 2 ) - 1 . « 2 ( 3 ) - ® 2 « = ^ jf = ^ = ^ = ^ ' '' 

The coefficients V are given in table H. 
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TABLE IV.- MATRICES * 9 {ffi,n g (P), ■ ■ ■ OCCURRING IN EQUATIONS (32) AND (33) 

Matrices ^g(j3),n fl (j3), . . . can be written in terms of nine submatrices, each 
of order p + 1, in the form 



©l03) 

© 2 (© 

© 3 (0 

:M**> = 

© 4 (/3) 

§5^ 

©6 {0) 


© 7 03 ) 

O g 03) 


General terms of matrices 

0 1 (^), © 2 (/3), . . . 

are given 


row and i to column. 


* g 03) 

m 

J g (P) 

Dz^ 

F Z W) 

T z (« 

.(IS) 


© t (/3) 

V(2,k,i) 

0 

0 

V(4,k,i) N6g(j3) 

V(l,k,i)n 2 - V(3,k,i)N2 
+ V(5,k,i) 

0 

V(3,k,i) 

0 

0 

V(4,k,i)m5g(j3) 

V(l,k,i) R 2 - V(2,k,i) m 2 
+ V(5,k,i) 


© 2 (|3) 

0 

0 

V(6,k,i) 

0 

0 

V(7,k,i) N6 6 ()3) 

0 

0 

0 

V(7,k,i) 

-V(6,k,i)m6g(0) 

0 


0 3 (/3) 

0 

0 

0 

V(8,k,i) 

-V(9,k,i)N6g(0) 

0 

0 

0 

V (9 ,k,i) 

0 

0 

V(8,k,i) m5g(/3) 


e 4 W) 

o 

o 

V(10,k,i) 

0 

0 

V(ll,k,i)N5 6 0) 

0 

0 

0 

V(ll,k,i) 
V(10,k,i) m6g(/3) 

0 


g 5 ^ 

V(12,k,i) 

0 

0 

V(14,k,i)N6g( J 3) 

V(l,k,i) R 2 - V(13,k,i) N 2 
+ V(15,k,i) 

0 

V(13,k,i) 

0 

0 

-V(14,k,i) m6 6 (|3) 

V(l,k,i) R 2 - V(12,k,i) m 2 
+ V(15,k,i) 


©g(/3) 

0 

V(16,k,i) 

-V(18,k,i) N6g(0 
0 
0 

-V(17,k,i)N 2 + V(19,k,i) 
0 

V(17,k,i) 

V(18,k,i) m6 6 (j3) 

0 

0 

-V(16,k,i) m 2 + V(19,k,i) 


e 7 0 3) 

o 

0 

0 

V(20,k,i) 

V (21 ,k,i) N5g(/3) 

0 

0 

0 

V(21,k,i) 

0 

0 

V(20,k,i) m8g(0) 


e 8 (ffl 

o 

V(16,k,i) 

V(18,k,i)N6g(/3) 

0 

0 

-V(17,k,i)N 2 + V(22,k,i) 
0 

V(17,k,i) 

-V(18,k,i) m6g(j3) 

0 

0 

-V(I6,k,i)m 2 + V(22,k,i) 


e 9 0S) 

V(23,k,i) 

0 

0 

-V(25,k,i)N6g(/3) 

V(l,k,i) R 2 - V(24,k,i) N 2 
+ V(26,k,i) 

0 

V(24,k,i) 

0 

0 

V(25,k,i)m6g((3) 

V(l,k,i) R 2 - V(23,k,i) m 2 
+ V(26,k,i) 


Notes: 8g(l) = 1; fig(2) = -1. 

The coefficients V are given in table II. 




TABLE V.- EQUATIONS FOR f 


fee(i), f 00( 2 )» 

f 00(l) = 


can be expressed in the form 


n r (l) - v r (2)-^ - 8Z 2 90 9Z 


UrO 


u 0 (1) ae + u e (2) 8Z 


'00 


( 1 )^ + u z^'S7; U z0 


90 


e r (l) 

u r (2) 

e r (3) 

e r (4) 

ue(l) 

t)0(2) 

u z (l) 

«z(2) 

B 22 ( 5 > 1 ) 

B 22 ( 4 ’ 2 ) 

B 2 3(8,2) 

B26( 17 > 2 ) 

B 22 ( 7 »l) 

B26(l°,l) 

626(8,1) 

623 ( 8 , 1 ) 1 

B 2 2(2,2) 

B 2 2(^ 3 ) 

623(8,8) 

626(18,3) 

B 2 2(4,2) 

B 2 6 ( 7 , 2 ) 

626(2,2) 

623(8,2) 

622(8,!) 

B 22 ( 2 > 2 ) 

623(6,2) 

B26( 13 >2) 

622(8,!) 

626(8,1) 

626(3,1) 

623(6, 1) 

623(6,!) 

B23(5,2) 

B33(9,2) 

1 Bgg(16,2) 

1 623(8,1) 

636(18,1) 

636(8,1) 

633 ( 8 , 1 ) 

B23( 6 . 2 ) 

623(8, 3 ) 

B33(9,3) 

B3g(16,!3) 

623(8,2) 

636(18,2) 

B 36 (6,2) 

633(9,2) 

B 2 6 ( 8 ,l) 

B26^ 7 ’ 2 ) 

B3g(15, 2 ) 

B 66 ( 14 > 2 ) 

B 26(l°'l) 

b 66(U,1) 

b 66^8,1) 

636(15,1) 

b 26( 3 ’ 1 ) 

B 26 (2, 2 ) 

B3e(6,2) 

666^ 3j ^ 

B26( 5 »l) 

666(8,1) 

b 66( 3 ,1) 

B36( 6 ,l) 

B26(13,2) 

j B26( 1 8,3) 

b 36(16,3) 

B 66 (12,3) 

B 26 (l 7 ’ 2 ) 

B 66 (l 4 , 2 ) 

666(13,2) 

636(16,2) 


fpr (1) =B P (9>1) t2 Ur0 


8t 

3 g2 

f - -B„(7 3) — U r o + Bp(10,2)— ^U 0O 

tprW - P V ’ 8t 2 90 P 9t 2 

2 

fpr(3) = - B P< 9 " 3 >^r^ Ur0 + B P (9>2) ^ Uz0 

tpe(D = - B P (10 ’ 2) iT^ Ur0 + B P (11 ’ 1) i2 Ueo 

fpzd) - - B P (9 ’ 2) ^JlS Ur0 + B P (9 ’ 1) I[2 Uz0 

The B terms are pure functions of the material properties and the inner and outer radii of the individual 
plies. 


- . .;B p (18,«) 

V pa(j+l) f 1 l 1. 1 . 1 1 . R ._- JL : R ; _SL; -3-s; 

- 4 4 U ) \i^ ; SO? 5’ la?’ «m’ ’ a a> 2 a(H 2 aa> 2 

[r + a(l)l 2 0 R a a(l3 Rfa + a(lj] . l£+XlL tA±X^\ c„(j) [R - ad)]'* 1 dR 

L ^^2 ; L R5(I)- ; ’ ' XD’ ad) ’ a(l)2 Rad ) 2 J 


p = 22, . . 66, p; when p = P, c /3 ~ P- 



pier 


TABLE VI.- FREQUENCY PARAMETER X FOR THREE-PLY CYLINDERS 
c ll(i): C 12 (i): Cl3(i) :C 22( i ) : C 23W : c 33( i ) : C 44 W : c 55( i ) : C 66< i ) 


= 0.08: 0.05: 0.07: 0.19: 0.32: 1: 0.04: 0.03: 0.34 (i = 1,2,3) 


_ C 33( 1 )/ C 33( 2 ) = C 3 3( 3 )/ C 33< 2 ) = 20 


(a) Thin cylinder; total thickness, 5 percent of outer radius; a(l) = 0.95; a(2) = 0.955; and a(3) - 0.995 


m 

First frequency 

Second frequency 

Third frequer 

icy 

Exact 

Thin 

shell 

Refined 

laminate 

Exact 

Thin 

shell 

Refined 

laminate 

Exact 

Thin 

shell 

Refined 

laminate 

0 

0.32461 

0.34393 

0.32475 

1.8186 

1.8323 

1.8321 

3.0037 

3.0582 

3.0531 

1 

.33631 

.36091 

.33646 

1.7031 

1.7141 

1.7139 

3.1438 

3.2082 

3.2024 

2 

.36737 

.40749 

.36753 

1.4689 

1.4757 

1.4756 

3.4615 

3.5521 

3.5440 

3 

.40447 

.47064 

.40467 

1.2612 

1.2657 

1.2654 

3.8451 

3.9760 

3.9636 

4 

.42507 

.52990 

.42533 

1.1822 

1.1875 

1.1857 

4.2502 

4.4368 

4.4173 

0 

.28282 

.28342 

.28286 

.91444 

.91614 

.91612 

1.5441 

1.5508 

1.5503 

1 

.30591 

.30722 

.30596 

.73851 

.73938 

.73936 

1.7807 

1.7917 

1.7909 

2 

.30838 

.31225 

.30846 

.60334 

.60419 

.60379 

2.2071 

2.2292 

2.2275 

3 

.21959 

.23314 

.21974 

.79860 

.80118 

.79965 

2.6748 

2.7153 

2.7122 

4 

.20414 

.24179 

.20444 

1.1413 

1.1467 

1.1445 

3.1559 

3.2241 

3.2187 

0 

.20999 

.20995 

.21000 

.22901 

.22904 

.22904 

.51742 

.51743 

.51757 

1 

.07054 

.07054 

.07054 

.41947 

.41954 

.41957 

.84022 

.84128 

.84126 

2 

.03594 

.03679 

.03602 

.78831 

.78949 

.78928 

1.3512 

1.3559 

1.3556 

3 

.06940 

.07419 

.06966 

1.1918 

1.1964 

1.1956 

1.8985 

1.9113 

1.9107 

4 

.12237 

.13560 

.12283 

1.5963 

1.6080 

1.6058 

2.4564 

2.4841 

2.4828 

0 

0 

0 

0 

0 

0 

0 

.40893 

.40875 

.40897 

1 

0 

0 

0 

.57734 

.57756 

.57764 

.59782 

.59819 

.59818 

2 

.02316 

.02375 

.02328 

.91106 

.91294 

.91266 

1.1935 

1.1964 

1.1963 

3 

.06372 

.06712 

.06401 

1.2851 

1.2910 

1.2900 

1.7847 

1.7946 

1.7943 

4 

.11781 

.12857 

.11829 

1.6695 

1.6833 

1.6808 

.2.3694 

2.3927 

2.3921 


(b) Thick cylinder; votal thickness, 20 percent of outer radius; a(l) = 0.8; a(2) = 0.82; a(3) = 0.98 


b 

a 

m 

First frequency 

Second frequency 

Third frequency 

Exact 

Thin 

shell 

Refined 
laminate; 
Pm = 1 

Exact 

Thin 

shell 

Refined laminate 

Exact 

Thin 

shell 

Refined laminate 

Pm = 1 

Pm = 2 

Pm = 1 

Pm ~ 2 

i 

0 

0.40438 

0.81033 

0.40466 

1.6205 

1.8399 

1.8312 

1.6637 

1.7475 

3.0636 

1.8508 

1.7901 


1 

.42140 

.90411 

.42173 

1.5294 

1.7059 

1.6938 

1.5648 

1.7530 

3.2348 

1.8543 

1.7947 


2 

.46495 

1.1402 

.46544 

1.3354 

1.4500 

1.4331 

1.3580 

1.7633 

3.6180 

1.8620 

1.803,7 


3 

.50904 

1.2187 

.50984 

1.1742 

1.4813 

1.2375 

1.1894 

1.7626 

4.0843 

1.8508 

1.8022 


4 

.52631 

1.1676 

.52774 

1.1540 

1.8484 

1.2171 

1.1695 

1.7309 

4.5893 

1.7982 

1.7670 

2 

0 

.31807 

.34461 

.31837 

.89129 

.91996 

.91833 

.89790 

1.4316 

1.5573 

1.5388 

1.4580 


1 

.35573 

.41939 

.35608 

.71061 

.72587 

.72295 

.71368 

1.4242 

1.8267 

1.4445 

1.4441 


2 

.33947 

.46341 

.33990 

.62782 

.68829 

.63563 

.62979 

1.3608 

2.3061 

1.3801 

1.3789 


3 

.28099 

.57619 

.28145 

.85051 

.98070 

.87356 

.85655 

1.5134 

2.8370 

1.5451 

1.5381 


4 

.33070 

.87921 

.33178 

1.1398 

1.3556 

1.2125 

1.1597 

1.7260 

3.3934 

1.7896 

1.7581 

8 

0 

.21844 

.21804 

.21861 

.22955 

.22999 

.22996 

.22965 

.54260 

.54266 

.54543 

.54332 


1 

.06638 

.06642 

.06639 

.46207 

.46637 

.46439 

.46267 

.86049 

.88568 

.88059 

.86560 


2 

.08773 

.13212 

.08787 

.82383 

.86418 

.84514 

.82944 

1.3288 

1.4482 

1.4065 

1.3499 


3 

.18459 

.33591 

.18496 

1.1547 

1.2943 

1.2350 

1.1765 

1.6937 

2.0573 

1.7601 

1.7266 


4 

.28616 

.61270 

.28705 

1.3692 

1.7296 

1.5441 

1.4177 

1.9282 

2.6843 

1.9793 

1.9503 

oo 

0 

0 

0 

0 

0 

0 

0 

0 

.44779 

.44467 

.44865 

.44802 


1 

0 

0 

0 

.60992 

.61894 

.61637 

.61160 

.64047 

.65076 

.64808 

.64238 


2 

.08004 

.11167 

.08019 

.91989 

.97362 

.95316 

.92872 

1.2280 

1.3015 

1.2853 

1.2418 


3 

.18123 

.31182 

.18159 

1.2068 

1.3758 

1.3065 

1.2338 

1.7339 

1.9523 

1.8120 

1.7733 


4 

.28393 

.58693 

.28481 

1.3884 

1.7937 

1.5749 

1.4405 

1.9938 

2.6030 

2.0327 

2.0090 





TABLE VII.- ACCURACY OF FREQUENCIES OBTAINED 
BY USING THIN-SHELL THEORY 


Thickness 
of laminate 

Nodal distance b 
Radius a 

Circumferential 
wave number, 
m 

Accuracy 


Thin 

Thick 

Long 

Moderate 

Short 

Reasonably 

accurate 

Inaccurate 

Highly 

inaccurate 

X 


X 

X 


^2 

X 



X 


X 

X 


>2 


X 


X 




X 

All 


X 



X 

X 



^1 

X 




X 

X 



>1 


X 



X 


X 

X 

All 



X 


Notes: 


1. Long: b/a > 8 

Moderate: 2 < b/a < 8 

Short: b/a < 2 

2. Reasonably accurate: error less than 5 percent 

Inaccurate: error 5 percent to 25 percent 

Highly inaccurate: error greater than 25 percent 

3. Thickness up to 5 percent of radius 
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TABLE VIH.- RELATIONSHIPS FOR d r , dg, AND d z AND FOR e Q , e 1( . . . 
(a) Relations for d r , dg, and d z 


j = 0 for — 

j > 0 for - 


1,2, 3, 4 

k = 5,6 

j is even and k = 1,2, 3, 4 or j is odd and k = 5,6 

j is even and k = 5,6 or 
j is odd and k=l,2,3,4 

1 

0 

( e 2 e 6 " e 4 e 5)/ e 7 

0 

e 0 

0 

( e 3 e 5 - e l e 6)/ e 7 

0 

0 

1 

0 

e 8/ e 9 

0 

0 

| e 2 e 6 + e 2 e 6 ' e 4 e 5 " e 4 e 5) e 7 ~ ( e 2 e 6 “ e 4 e 5) e 7]/ e 7 2 

0 

e 0 

0 

[( e 3 e 5 + e 3 e 5 ' e l e 6 “ e l e 6) e 7 ' ( e 3 e 5 ' e l e e) e 3 / e 7 2 

0 

0 

0 

0 

( e 8 e 9 " e 8 e 9)/ e 9 2 



( e 2 § 6 - e 4 g 5)/ e 7 

0 



( e 3®5 - e l S 6)/ e 7 

0 



0 

^8 /^9 



[( e 2 S 6 + e 2 g 6 ' e 4 g 5 _ e 4 e s) e 7 ' ( e 2 § 6 " e 4 § 5) e 7]/ e 7 2 

0 



[( e 3 g 5 + e 3 g 5 " e l g 6 ~ e l g 6) e 7 - ( e 3 g 5 ’ e l g 6) e 7]/ e 7 2 

0 


0 


(ege 9 - egegWeg 2 



— Continued 


TABLE VIII.- RELATIONSHIPS FOR d r , d fl> AND d z AND FOR e 0 , e v ■ ■ ■ 

(b) Relations for e 0 , ej, . - • 



Relation 


C n o(k ) 2 - C 55 m 2 - C 22 

e 0 

m(C 55 + C 2 2 ) " m(C 5 5 + C 12 ) «( k ) 

e l 

C n [a(k) + ]] 2 - C 55 m 2 - C 22 

e 2 

m|(C 55 + C 12 )[o:(k) + j] - C 55 - C 22 j 

e 3 

m ^-(C 55 + Ci 2 )[a(k) + jj - C 55 - C 22 ] 

e 4 

C 55 [a:<k) + j ] 2 ' C 22 m2 " C 55 

e 5 

(pfi 2 a 2 - C 44 N 2 ) d r (j -2, k) + Nd z (j-l,k) JJ(C 13 + C 44 )[a ! (k) + J - lj - C i3 + C 2 ^ 

e 6 

(pfi 2 a 2 - C 66 N 2 )d e (j-2, k) + mNd z (j-l, k)(C 23 + Cqq) 

e 7 

e l e 4 ' e 2 e 3 

e 8 

{{C 44 + C 13 ) [o(k) + j - 1] + C 44 + C 23 ^Nd r (j-l,k) + (C 66 + C 2 3 )mNd 0 (j-l,k) + (p« 2 a 2 - C 33 N 2 )d z (j-2, k) 

e 9 

^ 66 m ^ " < -'44G 3r ^ + l ] 2 — — — 


Note: d r (-l,k), d 0 (-l,k), and d z (-l,k), when they occur, are zero 


(c) Relations for e^, e^, . . . 



Relation 


c 55 + c 12 

e 0 

C 55 + c 22 

e i 

2Cn[a;( k ) + j] 

e 2 

m( c 55 + c 12) 

e 3 

-m(C 55 + C 12 ) 

e 4 

2C 55 [Q:(k) + j] 

e 5 

( p fi2 a 2 - C 44 N 2 )d^(j-2,k) - N(Ci 3 + C 44 )d z (j-1, k) + Nd^Q-l, k) ^(C i3 + C 44 )[“( k ) + * ' 3 ' Cl3 + C23 j 

e 6 

(p«2 a 2 - C 66 N2)d' e (j-2,k) + (C 23 + C 66 )mNd^j-l, k) 

e 7 

e l e 4 + e l e 4 ‘ e 2 e 3 ' e 2 e 3 

e 8 

( pSi 2 a 2 - C 33 N2)di( j -2,k) * (C 44 * C 13 )Nd r (]-l, k) + {(C 44 ♦ C 13 )[»(k) + j - lj + C 44 + C 2 ^ Nd^(j-l.k) + (C 66 * C 23 )n,Nd’ 9 U-l. « 

e 9 

-2C 44 [a!(k) + 3 ] 


Note: d^(-l,k), d^(-l,k), and d^(-l,k), when they occur, are zero 
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TABLE VTU.- RELATIONSHIPS FOR d r> dp, AND d z AND FOR eo, ei, ... - Concluded 


(d) Values for §5, eg, • ■ • 



3 = /3 (j is even) 

j = /3 + 1 (j is odd) 

j g /3 + 2 

e 8 

e 8 



e 9 

044(30(6) - a] at a = a(6); 
thus, eg = 20440; (6) 



S 8 

e 8 



®9 

-C44 



e 5 


Nd z (/3,6) |J-(C 13 + C 44 )[o(6) + j - lj - C 13 + C 23 j 


e 6 


mNd z (^,6) (C 2 3 + Cgg) 


e 5 

(p« 2 a 2 - C 44 N 2 )d r (/3-l, 6) - Nd z (0,6) (C 13 + C44) 

+ Ndz(p,6) ^-(C 13 + C44) (a (6) + 3 - l] - Cj 3 + C 23 j 


S 6 


(pn 2 a2 - C 66 N 2 )d 0 (/3-l,6) + mNd;(/3,6) (C 23 + Cgg) 


a 8 



(p£> 2 a 2 - 033^)3^0-2,6) + ^(C 4 4 + C 13 )[o(6) + 3 - l] 
+ C 44 + C 2 ^Nd r (j-l, 6) + (Cgg + C 23 ) mNdp(j -1,6) 

e 9 



e 9 

g 5 



(pfl 2 a 2 - C 44 N 2 )d r (j-2, 6) + N^(C 13 + C 44 )[a(6) + j - l] 
- C 43 + C 23 j d z (j -1, 6) 

g 6 



(pfi 2 a 2 - C 6 gN 2 )d 0 (j-2, 6) + mN(C 23 + Cgg)d z (j-1, 6) 

B 8 



(pn 2 a 2 - C 33 N 2 )d z (j-2, 6) + (C 44 + Ci3)Nd r (j-l, 6) 

+ j(C 44 + Ci 3 )(o(6) + j - lj + C44 + C 23 j Ndr(j-1, 6) 
+ (C66 + C 23 )mNdg(j-l, 6) 

®9 



eg = -2044(0(6) + jj 

S 5 



(p« 2 a 2 - C 44 N 2 )d(.(j-2, 6) - (Ci3 + C 44 )N<^(J-1,6) 

+ n|-(C 13 + C 44 )(q'( 6) + j - l] - C 13 + C 23 j d4(j-l,6) 

®6 



(pfi 2 a 2 - CggN 2 ) dg(j -2, 6) + mN(C 23 + Cgg)d z (;j-1, 6) 
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